Spectral function of the electron in a superconducting RVB state 
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We present a model calculation of the spectral function 
of an electron in a superconducting resonating valence bond 
(RVB) state. The RVB state, described by the phase-string 
mean field theory, is characterized by three important fea- 
' tures: (i) spin-charge separation, (ii) short range antifer- 
romagnetic correlations, and (iii) holon condensation. The 
results of our calculation are in good agreement with data 
' obtained from Angle Resolved Photoemission Spectroscopy 
t (ARPES) in superconducting Bi2Sr2CaCu208+6 at optimal 
doping concentration. 

The spectral function of an electron, A(k,a;), being 
the probability of finding an electron with momentum k 
and energy w, is a fundamental quantity in any descrip- 
tion of interacting electrons. ARPES is a very power- 
ful and direct experimental technique used to measure 
the spectral function in electronic systems. In the last 
decade, ARPES measurements have contributed much 
to our understanding of the high temperature supercon- 
ductors (HTSC). Measurements done in the normal and 
the superconducting state of the HTSC reveal a variety 
of interesting features 0|, For instance, the normal 
state spectra are extremely broad, with the broad peak 
evolving into a hump at Tc- In addition to the hump, 
a very sharp peak appears below Tc at low binding en- 
ergies. This is observed most clearly near the Brillouin 
zone boundary around the M points, (0, ±7r) and (±7r, 0) 
[^. The strength of the sharp peak appearing below 
is proportional to the superfluid density [Q, Along 
the direction (0,0) — > (tt, tt), there is some debate if the 
data below Tc can be interpreted as showing a clear break 
between a coherent quasiparticle part and the broad in- 
coherent background (D, fft]. 

It was first noted by Anderson that spin-charge sepa- 
ration may provide a natural explanation for the ARPES 
results 0. The basic idea is that the hole created by a 
photon decays into a spinless holon and a neutral spinon 
excitation. Such a decomposition can explain the broad 
incoherent background seen in photoemission spectra, as 
well as the absence of a sharp quasiparticle peak in the 
normal state. The emergence of a sharp quasiparticle 
state below Tc is attributed to the condensation of holons. 
However, as we shall demonstrate in this paper, the inclu- 
sion of short range antiferromagnetic (AF) correlations 
induced by the photoholc is crucial to the understanding 
of the ARPES results. This has been pointed out in the 
context of photoemission from a Mott insulator p|. In 



this Letter, we consider the photoemission spectra ob- 
served in the superconducting state of optimally doped 
Bi2Sr2CaCu208+5 (Bi 2212). We present a model calcu- 
lation of the spectral function in a superconducting RVB 
state, incorporating three important features: (i) short 
range AF correlations, (ii) spin-charge separation, and 
(iii) holon condensation. We show that such a descrip- 
tion provides a consistent explanation for many of the 
generic features observed in the superconducting state of 
optimally doped Bi 2212. 

AF oscillations: We begin with the decomposition of 
the electron operator in the so-called phase string for- 
mulation, Cic — h\aic, where h\ is the bosonic holon 
creation operator, and aia is a composite spinon opera- 
tor that satisfies fermionic anticommutation relations . 
Unlike the conventional slave-boson formalism, the (com- 
posite) spinon operator a is defined by Uia = fcio-e'®'", 
where hi^ is the bosonic elementary spinon annihilation 
operator and , a nonlocal phase string operator whose 
origin lies in the fact that a hole moving through a locally 
AF background always picks up a string of sequential ± 
signs Q . The factor Qi„ imposes anticommutation rules 
on fljo- as well as Ci„. 

In the superconducting phase, the holons {h\) con- 
dense. An elementary nodal (d-wave) fermionic excita- 
tion emerges in the spinon sector and its creation oper- 
ator is related to the a-operator, to leading order as 
Okcrl^'G) Of — crwk7kCTl^G)j where Vk is defined in the 
usual manner as Wk = sgn(Ak)(l — Ck/£'k)^^^/V2, with 
-^k — v^^k + ^k ^'^^ — Ao{cos k^a — cos kya). Thus, 
the superconducting ground state, |^'g): in the phase 
string description looks quite similar to the d-wave mean 
field state in the slave-boson approach |11|. However, 
the above considerations (obtained by using equations of 
motion), miss the local AF structure around a hole that 
is created. To see this, let us consider a bare hole state 
created by the c-operator, Ciil^c) and measure the spin 
configuration around the hole by evaluating the quantity 
{^g\4^S^c,^\^g)- Now, 

a 

= -l\(^G\bl/^,\^G)\' + i|(*G|6,V,T|*G)P . (1) 

Using the mean field solution jl^, we get 
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m 

(2) 

and 

(vI/g|6,V,tI*g) = [(-1)'-^' + 1] E«™<t(0^™tO') ■ 

m 

(3) 

In the above, Um and Vm are the "coherence factors" 
of the mean field theory determined self consistently for 
a given hole concentration 5, within a Bogoliubov-de 
Gennes scheme, and lOmtj, a one particle wave function. 
Equations ^ and (||) are used in (|l|) to determine the 
spin density around the hole site. The results are shown 
in Fig. 1, where we have plotted the spin density as a 
function of the hole distance along the x axis for hole 
concentrations (5 = 0, and (5 = 1/7 ~ 0.14. As seen in the 
figure, the hole is surrounded by AF oscillations in the 
local spin density. The oscillation is to be expected from 
(|l|) and is indicative of local AF correlations. The ground 
state |\1/g) is a singlet and the presence of Ciilvfc) (an 
up spin hole) is accompanied by a spinon excitation with 
= —1/2. Assuming the (up spin) hole to be located on 
an odd sublattice site, our results show that the (down) 
spinon is created only on the even sublattice sites. Note 
the presence of nonvanishing (up) spin density on the 
odd sublattice sites, which we interpret as overscreening. 
Thus, the combination of the up spin hole on an odd site 
and the down (up) spin response on the odd (even) sites 
represents a spin-polaron with zero net spin. 

To account for this effect in determining the spectral 
function, we adopt the following empirical approach. We 
write 

Oiil^c) — '^Vj{'i')o:]i\'^G) + higher order terms , (4) 

j 

where rjj (i) ^ only for i and j not belonging to the same 
sublattice, and Oko- = |wk|7kCT (the d-wave sign in Uk will 
be absorbed into rj). We only retain the nearest and 
third nearest neighboring sites in the expansion (^; viz., 
Vj{^) = +i-)Vo/2, ifj^i± x{y); r,j{i) = +(-)r/i/2, if 
j — i±2x{x)±y{2y); rjj (i) « for all other sites. The sign 
of r]j{i) is from the d-wave symmetry of spinon pairing. 
Transforming to momentum space, we get, within this 
approximation, aktl^c) ~ ?7k|wk|7lki I*g) , where 

77k = 770 (cos k^a — cos kyu) 

+ 2771 (cos fc^jO cos2kya — cos kyU cos2kxa) . (5) 

Here we take 771/779 ~ 0.3 for 6 — 0.14. 

We may now ask what happens when an electron is 
created (as, for example, in inverse photoemission), i.e., 
cL|5'g) . It can be shown in this case that the t spinon 



is created at the site i with some residual amplitude ex- 
tended over other sites of the same sublattice. Neglecting 
the residual amplitude, to leading order we get 

"ktl^G) -"k7ktl*G) , 

which is essentially the same as the slave-boson result. In 
the above, = 1 — v^. Therefore, the main distinction 
between the usual d-wave slave-boson mean field theory 
and this calculation is the momentum dependent factor 
77k that arises in the hole channel (corresponding to the 
creation of a photohole in ARPES). The remainder of 
this paper is devoted to a study of this channel. The 
differences between particle and hole spectroscopies as 
well as a detailed comparison with slave-boson theories 
will be presented elsewhere. 

Spin-charge separation and holon condensation: In 
the superconducting state, the spinous are paired with 
d-wave symmetry as discussed earlier. The effec- 
tive hamiltonian for the holons is given by Hh — 
—th e^*-^ h\hj -[-h.c., where represents the gauge 
field due to the spinous seen by the holons. Since the 
spinous are paired, the mean field solution leads to the 
result X]n ^/j ~ 10 ■ Thus the superconducting 
state is described by paired spinous and a Bose conden- 
sate of holons that experience a tt flux around an ele- 
mentary plaquette. Based on our earlier discussion, we 
obtain the spectral function of the spinous, pa(k, lj) = 
^k"k'^(^ + ^k) ■ The spectral function for the holons is 
easily derived as 

p„(k, oj) = cos2 ^-^5{oj - et ) + sin^ ^<5(c. - e^^) , 

where cos6'k — 7k/(V2Ak). Here, 7k = cos fc^jC -I- cos fc^a, 
Ak ~ \J cos^ fc^o -I- cos^ fcj^a, and ejj_|_ = ±2t/iAk — /i/i. At 
T = 0, the holon chemical potential [i^ = — 2t/iAk=o = 
-2V2th. 

We now construct the spectral function of the elec- 
tron by using the operator decomposition of the elec- 
tron, da- — hoaia- + c[g., where < hi >= ho, describes 
the Bose condensate of holons, and c^^ =: /i| : Qia, with 
: hi := hi — ho. The spectral function of the electron 
yl'i(k, u), is expressed as a convolution of the spinon and 
holon spectral functions. It is easy to see that 

k 

(6) 

where pa and pf^ are the spectral functions of the spinon 
and holon respectively. For obvious reasons, we write the 
spectral function as the sum of an incoherent and a co- 
herent part, A'L = A^_ + A'L- Clearly, these two terms 
correspond to the two terms in the electron decomposi- 
tion. 
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The coherent part of the spectral function is obtained 
from the contribution of the holon condensate. We get 

Ali^^.Lo) = plilvl5{u + E^) , (7) 

where p'j^ (x 5 denotes the density of the holon conden- 
sate. This contribution is dubbed "coherent" , since it 
is a sharp peak appearing below (the temperature at 
which the holons condense) . The incoherent part of the 
spectral function is the convolution, 

A'L (k, c.) = ^ E '^k'+k^^'+k[cos^ + El,^^ + 

k' 

+ sin2^<5(w + i?^,+k + ek'+)] ■ (8) 

The prime in the summation indicates that the contribu- 
tion from the holon condensate is removed. 

Comparison with ARPES results: We are now in a po- 
sition to plot (k, Lo) for various momenta and compare 
with results from ARPES. For the calculation of the spec- 
tral function, we choose the following parameters. The 
dispersion of the spinous, ^k, is determined by the near- 
est and next nearest neighbor hopping integrals, ti = 75 
meV, t2 — 20 meV. The chemical potential, /i = — 62 
meV is chosen to mimic the topography of the observed 
Fermi surface. We choose a value of A = 20 meV for 
the gap, and th — 4ii — 0.3 eV for the holon dispersion. 
Before proceeding to discuss the results for A'i_ (k, lo) , we 
note that the following can be anticipated: (a) The coher- 
ent part of the spectral function, (k, uj) (x (x 6. So, 
the strength of the coherent peak will be proportional to 
the superfluid density, as observed experimentally |^ , . 
(b) As is evident from our results, ^l(k, cj) is strongly 
momentum dependent, owing to the factor 77^. This con- 
tribution is strongest around the M point and weakest 
around the F point, (c) The incoherent part of the spec- 
tral function is expected to produce a broad background. 
Such a broad background is always seen in photoemis- 
sion. Since, in our calculation, the background arises 
from a convolution of spinon and holon spectral func- 
tions, we expect a rather weak momentum dependence 
of the broad background. This may have already been 
observed experimentally 

In Fig. 2, we show the results for A'L(k,oj) at a point 
on the Fermi surface, k — (0.4, tt). We see clearly that 
the total spectral response is the sum of the incoherent 
and coherent pieces of the spectral function. The sharp 
peak at lower binding energy corresponds to A'^ and the 
peak is located at an energy uj — E^. The broad feature 
seen in the figure is the contribution from A'. The inco- 
herent part exhibits a low energy edge around E^, and 
a broad peak ("hump") slightly above it. The origin of 
these features is very much the same as in the spectra of 
the undoped insulator In both cases, the low energy 
edge is determined by the dispersion of the spinon, E^. 
The origin of the hump lies in the local AF correlations 



embodied in the factor 77k- Let us consider (|[). The 
factor ?7k is maximum at the M point. As uj increases, 
k' 7^ terms in (H) have to contribute to the sum. How- 
ever, for such terms, the factor ?7k+k' decreases, thereby 
causing a hump in the spectral funtion. The position 
of the hump is shifted from the edge by ej^,„. A simi- 
lar effect occurs in the photoemission of a single hole in 
the Mott insulator, where a broad hump arises from the 
coherence factors of the Schwinger boson mean field the- 
ory. We emphasize that the factor 77k in the present case 
and the coherence factors in the Schwinger boson theory 
of the insulator reflect the antiferromagnetic correlations 
that play an important role in the formation of the ob- 
served humps. Seen in this perspective, the hump is not 
directly related to the 41 meV resonance observed in the 
neutron scattering fl^ ], as conjectured in the literature 
10 though AF correlations clearly play a crucial role in 
both explanations. 

In Fig. 3, we show the behavior of A°_ (k, uj) for vari- 
ous k points. In the left panel, we show how the spectral 
function evolves as one moves away from the AI point in 
two perpendicular directions. As can be seen in the fig- 
ure, the peak-dip-hump structure gets more pronounced 
near the M point. In the right panel, we show the spec- 
tral function for a series of points on the Fermi surface. 
Again, we see that the factor ryk causes the coherent peak 
as well as the hump to diminish away from the M point. 
In particular, along the direction F — > (tt, tt), these fea- 
tures are absent in our calculation. Though this is an 
approximate result, it shows that the inclusion of AF cor- 
relations suppresses the peak-dip-hump structure along 
this direction. 

To conclude, our model calculation of the spectral 
function of an electron in a superconducting RVB state 
explains many of the generic features observed in the 
photoemission spectra of optimally doped Bi 2212. We 
find that the inclusion of short range antiferromagnetic 
correlations is extremely crucial to the understanding of 
these features. In particular, we show that the presence 
of a hole induces Fricdel-likc oscillations in the local spin 
density surrounding the hole. We propose a scheme to 
incorporate this effect in the calculation of the spectral 
function. When this is taken into account in conjunc- 
tion with spin-charge separation and holon condensation 
below superconducting T^, many of the features seen in 
photoemission can be explained naturally. Our calcula- 
tion shows that these features have analogues in the pho- 
toemission spectra of the undoped Mott insulator where 
the short range antiferromagnetic correlations play a very 
important role too. 
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Fig. 1 Spin configuration near an up-spin hole, created 
by Ci on the ground state at 5 = and 0.14, along the x 
axis. 



Fig. 2 The spectral function at k =(0.4, tt), which is 
a Fermi point at the Brillouin zone boundary (see the 
inset). The dashed curve shows only the incoherent part 
whose low energy edge coincides with the coherent peak 
position while its 'hump' lies slightly above the edge. 

Fig. 3 The spectral function at different momenta 
marked by full circles in the insets of two panels. The 
left panel corresponds to scans near the M point along 
two directions, while the right represents momenta on 
the Fermi surface. 
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